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Abstract. 

We review results from general relativistic axisymmetric magnetohydrody- 
namic simulations of accretion in Sgr A*. We use general relativistic radiative 
transfer methods and to produce a broad band (from millimeter to gamma-rays) 
spectrum. Using a ray tracing scheme we also model images of Sgr A* and com- 
pare the size of image to the VLBI observations at 230 GHz. We perform a 
parameter survey and study radiative properties of the flow models for various 
black hole spins, ion to electron temperature ratios, and inclinations. We scale 
our models to reconstruct the flux and the spectral slope around 230 GHz. The 
combination of Monte Carlo spectral energy distribution calculations and 230 
GHz image modeling constrains the parameter space of the numerical models. 
Our models suggest rather high black hole spin (a* ~ 0.9), electron temperatures 
close to the ion temperature (Ti/T e ~ 3) and high inclination angles (i ps 90°). 



1. Introduction 

Observations of the Galactic Center provide strong evidence for the existence 
of a supermassive black hole in Sgr A* (which hereafter refers to the black 
hole, the accretion flow, and the radio source). Sgr A*'s proximity allows us to 
perform observations with higher angular resolution than other galactic nuclei. 
Estimates of Sgr A*'s mass M = 4 5±0. 4x 10 6 M m and distance D = 8.4±0.4 kpc 
( Ghez et "ail 120081 . iGillessen et al.l 120091 ) indicate that it has the largest angular 



size of any known black hole (GM/(c 2 D) ~ 5.3/zas). 

Recent 230 GHz VLBI constrains the struct ure of Sgr A* on angul ar scales 
comparable to the size of the black hole horizon ( Doeleman et al. 20081 . see also 



Doeleman contribution to this conference proceeding). Using a two-parameter, 
symmetric Gaussian brightness distribution model VLBI infers a full width at 
half maximum FWHM = 37l}g //as. This is smaller than the apparent diameter 

of the black hole: ~ 2\^27GM/(c 2 D) ~ 55/ias (this depends only weakly on 
black hole spin). Sgr A* radio - submm emission is usually modeled as syn- 
chrotron emission, with the turnover at ~ 230 GHz indicating transition from 
optically thick to optically thin emission. The model of accretion and its geom- 
etry is still under debate because Sgr A* is dim at shorter wavelengths (in NIR 
and X-ray band) or completely obscured (UV and optical light). Moreover in 
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the NIR and X-ray the source is not resolved and we only have upper limits for 
its quiescent luminosity. 

Sgr A* is a strongly sub-Eddington source (LboL ~ W 9 LEdd)- Models sug- 
gest that the source is accreting inefficiently (in the sense that L bol /(Mc 2 ) < 1), 
which justifies treating the dynamics of the flow and its radiative properties 
independently. Therefore most of Sgr A* models published to this day con- 
sists of two separate parts: plasma dynamics model and/or radiative transfer 
model. We can further categorize the dynamical models into: accretion flow 
vs. outflow models, stationary vs. time-dependent, Newtonian/post-Newtonian 
vs. General Relativistic, models covering small region (a couple of gravitational 
radii) vs. large region (thousands of gravitational radii). The radiative transfer 
modeling usually uses ray tracing (always relativistic) or Monte Carlo methods 
(nonrelativistic as well as fully relativistic). Ray tracing allows one to model 
source images at sub-mm frequencies and the SED of direct synchrotron emis- 
sion (radio, sub-mm). Monte Carlo allows one to model the Compton scattering 
and multiwavelength (from radio to gamma-rays) SED of Sgr A*. In Table [T] we 
summarize the recent progress in models. 



Reference 


dynamical 


radiative 


plasma 


range 




model 


model 




of model 


Naravan et al. (1998) 


stat. rel. ADAF 


non-rel. MC 


th 


10°i?g 


Markoff et al. (2001) 


Jet 


scaling 


non-th 




Yuan et al. (2003) 


stat non-rel. RIAF 


non-rel rays 


th+non-th 


2 x 10 5 i? s 


Ohsuea et al. (2005) 


MHD-time dep. 


non-rel. MC 


th 


60Rg 


Goldston et al. (2005) 


MHD-time dep. 


polarized non-rel. rays 


th+non-th 


512R g 


Broderick & Loeb ( 2006b) 


stat. non-rel RIAF 


polarized RT 


non-th 


2 x W 5 R g 


Moscibrodzka et al. (2007) 


MHD-time dep. 


non-rel. MC 


th+non-th 


2.4 x 10 3 R g 


Loeb & Waxman (2007) 


Jet 


scaling 


th+non-th 




Huang et al. (2007) 


stat. RIAF 


RT 


th 


2 x lO 5 ^ 


Markoff et al. (2007) 


Jet 


non-rel rays /w corr 


non-th 




Huang et al. (2009) 


stat. rel. RIAF 


RT 


th 


lO 4 ^ 


Broderick et al. (2009) 


stat. rel. RIAF 


RT 


th+non-th 


2 x 10 r °R g 


Chan et al. (2009) 


MHD-time dep. 


non-rel rays /w corr. 


th+non-th 


AZRg 


Yuan et al. (2009) 


stat. rel. RIAF 


RT 


th 


IQORg 


Hilburn et al. (2009) 


GRMHD-time dep. 


non-rel MC 


th 


40R g 


Dexter et al. (2009) 


GRMHD-time dep. 


RT 


th 


AORg 


Moscibrodzka et al. (2009) 


GRMHD-time dep. 


RT + rel. MC 


th 


40R g 



Table 1. Summary of selected models of Sgr A*. Abbreviations: RT- 
ray tracing, MC-Monte Carlo, GR-general relativistic, RIAF-radiatively in- 
efficient accretion flow, ADAF-advection dominate accretion flow, plasma- 
particles distribution, th-thermal, non-th-non-thermal, range-model radial 
range. 



2. Methodology 



Our numerical model of accretion onto Sgr A* consist of: a physical model 
of the accretion flow dynamics with its numerical realization, and a radiative 
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transfer model. We assume that the accreting plasma is geometrically thick, 
optically thin, turbulent and it accretes onto a spinning black hole. The spin 
angular momentum J of the black hole, whose magnitude is parameterized by 
a* = Jc/(GM 2 ), is assumed to be aligned with the angular momentum of the 
accretion flow. 

The accretion flow in Sgr A* is collisionless. We assume that ions and elec- 
trons have thermal distributions, possibly with different temperatures. In our 
model we allow the electron temperature T e to differ from the ion temperature Tj, 
but we fix the ratio (Tj/T e = const). The plasma equation of state is described 
by the P = e(7 — 1) relation with 7 = 13/9 (non-relativistic ions and relativistic 
electrons). A more physical model would evolve T e and Tj independently with a 
model for dissipation and energy exchange between the electrons and ions, but 
this would complicate the model and introduce a host of new parameters. 

We rea lize the physical mod el using an axisymmetric version of the GRMHD 
code harm ( Gammie et al.l [20031 ). As initial conditions w e adopt an analytical 
mode l of a thick disk (a torus) in hydrostatic equilibrium ([Fishbone fc Moncriei 
1976). Since the code is not designed to evolve MHD equations in vacuum, we 



surround the equilibrium torus with a hot, low density plasma which does not 
influence the torus evolution. We seed the torus with a poloidal, concentric loops 
of weak magnetic field. Small perturbations are added to the internal energy 
which allows to for magnetorotational instability and turbulence development. 
We solve the GRMHD evolution equations of evolution until a quasi-equilibrium 
accretion flow is established (meaning that the flow is not evolving on the dy- 
namical timescale). The details of torus initial s etup and discussion of the flow 
evolution is presented in McKinnev &: Gammid (|2004l ). 

Our numerical domain extends from the black hole event horizon to 40GM/c 2 
1.8 AU or 210/uas, but is only in equilibrium to ~ IhGM/c 2 = 0.7AU or 80/ias. 
Since low frequency emission from Sgr A* is believed to arise at larger radius, 
we are unable to model the low frequency (radio, mm) portion of spectral en- 
ergy distribution (SED). Our (untested) hypothesis is that at r < YbGMjc 2 
the model accurately represent the inner portions of a relaxed accretion flow 
extending over many decades in radius. 

Our fluid dynamical model is scale free but the radiative transfer calculation 
is not. Therefore we need to specify the simulation length unit C = GM/c 2 , time 
unit T = GMj c 3 , and mass unit M. (which scales the mass accretion rate) . C 
and 7~ are defined by the black hole mass adopted from observations. M is a free 
param eter set to reproduce the submillimeter flux F230GHZ = 3.4Jy, (jMarrond 
2006). To "observe" the numerical model we need to specify the inclination % of 



the black hole spin to the line of sight. The distance to the observer is assumed 
to be D = 8Akpc. 

The SED is calculated in a Monte Ca rlo fashion. We use grmonty - a gen- 
eral relativistic radiative transfer scheme ( Dolence et al. 20091 ) to calculate the 
light propagation in the strong gravity taking into account synchrotron emis- 
sion and absorption and Compton scattering, grmonty has been extensively 
tested. As one of our tests we compared the code performance in a flat space 
with independent Monte Carlo scheme Sphere (kindly provided by its author, 
Shane Davies). In FigfTJwe present a standard test for Compton scatte ring in 
the spherical, homogenous, cloud of hot plasma (jPozdniakov et aT1ll979l ). 
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Figure 1. Code testing: grmonty (points) and Sphere (independent Monte 
Carlo code for radiative transfer, solid line) SED comparison. In this test 
problem a thermal point source emits photons from the center of the spherical, 
hot plasma cloud. The cloud is characterized by two parameters: optical 
thickness r (here r = 0.1) and electron temperature O e (0 e = kT e /m e c 2 = 
4.). Lower panel shows the difference between two SEDs. 



In the present models we compute spectra through a time slice of a simula- 
tion as if the light had infinite speed ('fast light' approximation). In reality the 
light crossing time is comparable to the dynamical time which suggest radiative 



transfer through the changing medium (t^. 



20s in Sgr A*). Finally 



we use a ray-tracing scheme, ibothros (jNoble et al.l 12007). which accounts only 
for synchrotron emission and absorption in order to calculate the 230 GHz in- 
tensity maps of an accretion flow as seen by a distant observer (again in the fast 
light approximation). To estimate the size of the emitting region we calculate 
the eigenvalues of the matrix formed by taking the second angular moments of 
the image one the sky (the lengths of the principal axes). The eigenvalues along 
the major and minor axis are then compared to the FWHM VLBI constraint. 

We accept or reject numerical models based on several observational con- 
straints: the flux at 230 GHz (we scale our models to reconstruct the flux F=3.4 
Jy at 230 GHz), submillimeter spectral slope a around a tu rn over frequency 
(230-690 GHz, a changing from - 0.46 to 0.08 iMarrond I2006T) . upper limit for 
near-IR quiescent luminosity (e.g. iGenzel et al 20031 ). an upper limit for X-ray 
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luminosity L x < 2.4 x 10 33 [ergs ' 1 ] (|Baganoff et alJl2003l >. size of the image at 
230 GHz (jDoeleman et a l. 2008). We compare only time-averaged SEDs and 
time- averaged images to the observations. 



3. Results 
Best-bet model 

We vary three (free) model parameters: spin of the black hole (a* =0.5, 0.75, 
0.84, 0.94, 0.96, 0.98), the ion to electron temperature ratio (Ti/T e = 1,3,10), 



and inclination of the observer with respect to the spin axis (i = 5 ,45 , and 

o 

85 ). In Table [2] we show a set of accretion flow models in which SEDs are 
consistent with multiwavelength data as observed at at least one inclinations. 
We find the "best-bet" model, that best satisfies all the observational constraints 
(including VLBI 230 GHz size constrains) has spin a* = 0.94, Tj/T e = 3 and 

o 

i = 85 . All other models are inconsistent with the observed SED or only 
marginally consistent with the VLBI size measurement. 

In Figure [2] we show the SED of our best-bet model. The SED peaks 
around 690 GHz due to thermal synchrotron emission. Below 100 GHz it fails 
to fit the data because that emission is produced outside the simulation volume. 
A second peak in the far UV is due to the first Compton scattering order, at 
higher energies the photons are produced as an effect of two or more scatterings. 

In Figure [3] we show the dynamical model corresponding to our best-bet 
model. We show maps of number density density, magnetic field strength and 
electron temperature (6> e = kt,T e /m e c 2 ) on a single time slice. 

Figured] maps the points of origin for photons below and in the synchrotron 
peak (100-690 GHz), in the NIR (10 13 - 10 14 Hz), and in the X-ray (2-8 keV). 
The figure again corresponds to a single time slice from our best-bet model 
presented in Figures (2] [3] Most of the submillimeter emission originates near 
the mid plane at 4 < rc 2 /(GM) < 6. NIR photons are produced in the hot 
regions close to the innermost circular orbit (risco ~ 2.04GM/c 2 ). Photons in 
the X-ray band are produces mainly by scatterings in the hottest parts of the 
disk also close to the rjsco- A small fraction of photons are emitted from the 
funnel wall at large radii (15 — 40GM/c 2 ). 

Summary of parameter survey 

We can draw the following general conclusions from our work: (1) Very few 
of the time averaged SEDs based on a single-temperature (Tj/T e = 1) models 
produce the correct spectral slope between 230-690 GHz. The exception is edge- 

o 

on tori (i = 85 ) around a fast spinning black hole (a* = 0.98,0.96), but these 
models overproduce NIR and X-ray flux. 

o 

(2) For Ti/T e = 3 the only one model with a* = 0.94 at i = 85 is consistent 
with all observational constraints. For i = 85°, models with spins below a* = 
0.94 are ruled out by the inconsistent spectral slope, and models with higher 
spins (a* > 0.94), although consistent with the observed a, overproduce the 
quiescent NIR and X-ray emission. All models with Tj/T e = 3 observed at 

o o 

i = 5 and 45 are ruled out by the inconsistent a. 
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run 


a* 


i 

[deg] 


< M ■ 10~ 9 > 
[Mo/yr- 1 ] 


a 


lo-gio L x 

[ergs' 1 ] 


V 


cons, 
w/ obs.? 






5 


1.9 


-1.68 


31.5 


3.5 x 10" 2 


NO 




n 94 


45 


1 7 


-1 97 

L.Li \ 


31 8 


3.1 x 10~ 2 


NO 






85 


1.86 


-0.44 


32.9 


3.4 x 10~ 2 


YES 






5 


90.8 


-1.37 


30.1 


5.4 x 10~ 4 


NO 


AlO 


0.5 


45 


117.1 


-0.2 


31.4 


6.7 x 10~ 4 


YES 






85 


369.0 


1.38 


33.5 


1.6 x 10" 3 


NO 






5 


38.3 


-1.05 


30.0 


1.6 x 10~ 3 


NO 


BIO 


75 


45 


50 9 


04 


31 8 


2.0 x 10 -3 


YES 






85 


190.6 


1.49 


34.6 


6.9 x 10~ 3 


NO 






5 


19.5 


-1.15 


31.4 


5.1 x 10" 3 


NO 


ca o 


87^ 


45 


93 6 


-0 07 


3? 


6.2 x 10~ 3 


YES 






85 


41.4 


1.19 


34.2 


1.7 x 10~ 2 


NO 






5 


13.7 


-0.93 


31.8 


1.1 x 10" 2 


NO 


D10 


0.94 


45 


15.2 


-0.05 


32.3 


1.1 x 10~ 2 


YES 






85 


31.2 


1.17 


34.4 


2.5 x 10" 2 


NO 






5 


13.6 


-0.40 


32.5 


2.6 x 10~ 2 


YES 


E10 


0.97 


45 


14.3 


0.2 


33.1 


2.8 x 10~ 2 


NO 






85 


26.5 


1.19 


35.2 


5.1 x 10" 2 


NO 






5 


9.07 


-0.6 


32.7 


5.3 x 10~ 2 


YES 


F10 


0.98 


45 


9.16 


-0.08 


33.2 


5.2 x 10" 2 


NO 






85 


15.5 


1.12 


35.4 


9.0 x 10~ 2 


NO 



Table 2. Chosen set of models with temperature ratio 3 and 10 which 
SEDs are consistent with observations at least at one observation angle. The 
columns from left to right are: run ID (the number is the run name indicates 
Ti/T e ), dimensionlcss spin of the black hole, inclination angle of the observer 
with respect to the black hole spin axis, averaged rest mass accretion rate, a 
spectral slope between 230-690 GHz (F <~ z/*), and luminosity in the X-rays 
(at v <~ 10 18 Hz), the radiative efficiency 77 = Lbol/ Mc 2 , and whether the 
model is consistent with the data. We do not show models with Ti/T e = 1 
here because they are all inconsistent with the observations. 



o 

(3) For Ti/T e = 10, we find that all models with i = 85 are ruled out by 
both incorrect a and violation of NIR and X-ray limits. For lower inclination 

angles (i = 5°, 45°) a few models (E10 and F10 with i = 5°, and AlO, B10, C10, 
and D10 at i = 45°) reproduce the observed a. These models are consistent 
with X-rays and NIR limitations. Models E and F for i = 45° are ruled out by 
NIR and X-ray limitations whereas models AlO, B10, C10 and D10 for i = 5° 
produce a which is too small. 

(4) The dependence on a* arises largely because as a* increases the inner 
edge of the disk — the ISCO - reaches deeper into the gravitational potential of 
the black hole, where the temperature and magnetic field strength are higher. 
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Figure 2. Left panel: SEDs computed based on a single time slice t = 
1680GM/c 3 (thick line) (see Figure [3] for the distributions of the physical 
variables corresponding to the same time) along with the time averaged spec- 
trum (thin line) of our best-bet model. Observational points ar e taken from: 
Falcke et al]|1998l lAn et al.|[2005l iMarrone et aI1l2006l at r ad io. iGenzel et al.l 



2003| at NIR (1.65, 2.16, and 3.76 ^m) and iBaganoff et al.l 120031 at X-rays 
(2-8 keV) . The upper limit s in the NIR band a re taken fro m lMelia fc Falcke 



2001 



(30, 24.5 and 8.6 /mib lSchodel et al.ll2007l (8.6 fim) and lHornstein et al 



2007 (2 fim). The points in the NIR at flaring state are fromlGenzel et a! 



20031 (1.65, 2.16, and 3.76 H , and iDodds-Eden et al.l l2 009l C3.8 urn). An 
exam ple of X-ray flare (Lx = 1 x 10 35 ergs -1 ) is taken fromlBaganoff et al.l 
2001. Right panel: corresponding to SED in the left panel Image of the accre- 
tion flow at 230GHz correspindign to t — 1680GAf/c 3 . Intensities are given 
in units of er(7s _1 pixelsize _2 Hz _1 sr _1 , where the pixel size is 0.82/Ltas. The 
image shows inner 40GM/c 2 . The white circle marks FWHM=37 /ias of a 
symmetric Gaussian brightness profile centered at the image centroid. 



In the disk mid-plane, the temperature is proportional to the virial temperature 
and scales with radius e oc 1/r. B oc 1/r, while the density ~ r, below 
the pressure maximum (longer simulations in 3D show a more relaxed, declining 
density profile). Holding all else constant this implies a higher peak frequency for 
synchrotron emission, a constant Thomson depth (in our models the Thomson 
depth at the ISCO is roughly constant, since the path length 1/ ~ rjsco but 
the density ~ rjsco), and a larger energy boost per scattering A 166g, as 
can be seen in comparing models with different spin in Figure 4. The X-ray flux 
therefore increases with a* because e at the ISCO increases. 

(5) The dependence on Tj/T e is mainly due to synchrotron self-absorption, 

o 

which is strongest at high inclination. For example, because the i = 85 , 
Tj/T e = 10 model is optically thick at 230GHz the emission is produced in a syn- 
chrotron photosphere well outside rjsco- The typical radius of the synchrotron 
photosphere ranges between 15 GM/c 2 for low spin models (a* = 0.5,0.75) and 
8 GM/c 2 for high spin models (a* > 0.75). The 230GHz flux can then be pro- 
duced only with large A4; as A4 increases the optically thin flux in the NIR in- 
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Figure 3. Accretion flow structure in our best-bet model with a* = 0.94 
and with Tj/T e = 3 (model D3). The number density overplotted with the 
velocity field, the magnetic field strength, and the electron temperature are 
shown in the left, middle, and right panel respectively. The axis units are 
GM/c 2 . The figure shows a single time slice. 



creases due to increasing density and field strength. The scattered spectrum also 
depends on Ti/T e since the energy boost per scattering is ~ 160g oc l/(Tj/T e ) 2 . 

(6) The inclination dependence is a relativistic effect. M. is nearly indepen- 
dent of i (it varies by ~ 10%, except for Ti/T e = 10, which due to optical depth 
effects has much larger variation), so models with different inclination are nearly 
identical. Nevertheless the X-ray flux varies dramatically with i, increasing by 

o 

almost 2 orders of magnitude from i = 5° to i = 85 . This occurs because 
Compton scattered photons are beamed forward parallel to the orbital motion 
of the disk gas. The variation of mm flux with i is due to self-absorption. The 
mm flux reflects the temperature and size of the synchrotron photosphere. At 
lower i the visible synchrotron photosphere is hotter than at high i. 



4. Future prospects 

We have made a number of approximations that will be removed in future mod- 
els of Sgr A*. We plan to (1) add cooling, (2) run 3D (rather than axisym- 
metric) models, (3) model a larger range of radii, (4) include a population of 
nonthermal electrons, (5) eliminate the fast-light approximation by doing fully 
time-dependent radiative transfer. 
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Figure 4. Maps showing the point of origin for photons in our best-bet 
model. We show the logarithm of the sum of the photon weights in each zone, 
which is proportional to the number of photons seen at 100 — 690GHz (left), 
10 13 — 10 14 Hz (middle) and 2-8 keV (right) band. The axis scale units are 
GM/c 2 . The figure presents a single time slice at t = 1680GM/c 3 . Note that 
gray scale bands differ in scales. 
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